arXiv: 1508.04444vl [hep-ph] 18 Aug 2015 


ACFI-T15-11 


TeV Lepton Number Violation: From Neutrinoless Double /3-Decay to the LHC 

Tao Peng,^’[^ Michael J. Ramsey-Musolf,^’and Peter Winslow^’j^ 

^Department of Physics, University of Wisconsin-Madison, Madison, WI 53706 
^Amherst Center for Fundamental Interactions, Department of Physics, 

University of Massachusetts-Amherst, Amherst, MA 01003 USA 
^Kellogg Radiation Laboratory, California Institute of Technology, Pasadena, CA 91125, USA 

We analyze the sensitivity of next-generation tonne-scale neutrinoless double /3-decay {QvPfi) ex¬ 
periments and searches for same-sign di-electrons plus jets at the Large Hadron Collider to TeV scale 
lepton number violating interactions. Taking into account previously unaccounted for physics and 
detector backgrounds at the LHC, renormalization group evolution, and long-range contributions 
to Oirfip nuclear matrix elements, we find that the reach of tonne-scale QvjSP generally exceeds that 
of the LHC. However, for a range of heavy particle masses near the TeV scale, the high luminosity 
LHC and tonne-scale Oiz/3/3 may provide complementary probes. 


Total lepton number (L) is a conserved quantum num¬ 
ber at the classical level in the Standard Model (SM) of 
particle physics, yet it is not conserved in many scenarios 
for physics beyond the Standard Model (BSM). Among 
the most widely considered classes of BSM scenarios ad¬ 
mitting lepton number violation (LNV) is the see-saw 
mechanism [TH^ for neutrino mass generation that pos¬ 
tulates the existence of heavy right-handed Majorana 
neutrinos. Nr. The associated mass operator violates 
L by two units, and the coupling of the Nr with the 
left-handed SM lepton-doublet induces a corresponding 
L-violating Majorana mass for the three lightest neu¬ 
trino mass eigenstates. Experimentally, the observation 
of neutrinoless double beta-decay (0i//3/3-decay) of atomic 
nuclei would provide direct evidence for LNV. This ob¬ 
servation would also indicate the existence of a Majorana 
mass term for the lightest neutrinos[ 6 ], consistent with 
the prediction of the see-saw mechanism. 

Recent results from the EXO[7], GERDA[ 8 ], and 
KamLand- ZEN [i HU] experiments have placed stringent 
upper limits on the Ojz/3/3-decay half lives of ^®Ge 

and ^^®Xe on the order of a few times 10^® years. When 
interpreted in terms of the exchange of light Majorana 
neutrinos, these limits imply an upper bound of order 
100-400 meV on the Ojz/3/3-decay effective mass TO/ 3 / 3 , de¬ 
pending on the value of the nuclear matrix element em¬ 
ployed in this extraction[3S]. The next generation of 
“tonne scale” Ozz,5/3-decay searches aim for half life sensi¬ 
tivities of order > 10 ^^ years, with a corresponding TO/ 3 ^ 
sensitivity on the order of tens of meV, consistent with 
expectations based on the inverted hierarchy (IH) for the 
light neutrino mass spectrum. In this interpretive frame¬ 
work, a null result would imply that either neutrinos are 
Majorana particles with a mass spectrum in the normal 
hierarchy (NH) or that they are Dirac fermions. 

It is possible that neutrino oscillation studies may de¬ 
termine the neutrino mass hierarchy before the next gen¬ 
eration Ojz/3/3-decay searches reach their goal sensitivity. 
Should the hierarchy turn out to be normal, for which 
TO/ 3 ^ may range anywhere from the present upper bound 


to a magnitude <<10 meV, a null result from the tonne- 
scale Oiz/3/3-decay experiments would not be surprising. 
However, alternate decay mechanisms could still lead to 
observation of a signal in the next generation searches, 
even if the light neutrino spectrum follows the NH and 
the value of mpp is experimentally inaccessible. These 
mechanisms include radiative neutrino mass scenarios |llj 
and the low-scale see-saw mechanism [TUHHj [ISj • In these 
scenarios, the LNV interactions may involve particles 
whose masses are of order one TeV and whose exchange 
generates short range interactions that lead to the Ojz/ 3/3- 
decay. Straightforward arguments indicate that the re¬ 
sulting Ojz/3/3-decay half-life can be of order 10^^ yr or 
shorter, comparable to expectations based on the three 
light Majorana exchange mechanism and the IH[T7j. The 
associated light Majorana masses may nevertheless follow 
the NH with mpp well below the meV scale. 

How might one experimentally distinguish the TeV 
LNV scenario for 0zz/3,d-decay from the more conventional 
paradigm based solely on the exchange of light Majorana 
neutrinos? One possibility is to analyze experiments that 
search for charged lepton flavor violation, as discussed in 
Ref. |17] . Another, perhaps more direct, means is to 
search for the LNV interactions in high energy collider 
experiments (see, e.g. [I8H28] L 

This possibility has recently been explored by the au¬ 
thors of Refs. [UilUU], who utilized a simplified model 
framework to analyze the relative sensitivities of tonne- 
scale Ozz,5/3-decay experiments and searches for LNV sig¬ 
nals at the GERN Large Hadron Gollider. These authors 
performed a systematic classification of simplified mod¬ 
els that one may map onto more complete theories, such 
as R-parity violating supersymmetry. They find that in 
a broad range of cases the LHG with 300 fb“^ of inte¬ 
grated luminosity (corresponding to the end of Run H) 
would achieve substantially greater reach for TeV-scale 
LNV interactions than would the tonne-scale Ozz/3/3-decay 
searches [UU]. If verified, the prospective LHG exclusion 
of TeV scale LNV, coupled with identification of the 
NH for the light neutrino mass spectrum, could render 
the prospects for discovery with tonne-scale Ozz/3/3-decay 
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W-W--b2j 

W-Z-t2j 

ZZ-t2j 

W+2j 

tt 

tt 

t+3j 

W-+3j 

4j 


Before Cuts 

0.142 

0.541 

6.682 

0.628 

903.16 

68.2 

6.7 

0.45 

15.09 

362.352 

0.0038 

Signal Selection 

0.091 

0.358 

4.66 

0.435 

721.7 

28.9 

2.37 

0.22 

11.73 

72.03 

0.0031 

HrOets) > 650 GeV 

0.054 

0.04 

0.187 

0.015 

5.6 

0.266 

0.025 

0.0003 

0.102 

0.027 

0.0213 

> 130 GeV 

0.039 

0.029 

0.105 

0.008 

0.163 

0.127 

0.024 

3x10“'* 

0.101 

0.027 

0.0493 

< 40 GeV 

0.036 

0.005 

0.036 

0.007 

0.126 

0.014 

0.005 

3x10"® 

0.03 

0.017 

0.0684 

< 2.2 

0.033 

0.003 

0.022 

0.005 

0.093 

0.009 

0.004 

2x10“® 

0.019 

0.011 

0.0738 


TABLE I: Cut-flow designed for optimizing signal relative to backgronnd. Note: kinematic cuts are not commutative. 


searches less compelling than presently considered. 

In what follows, we revisit the analysis of Refs. [2^1^ 
and find that their conclusions regarding the LHC reach 
may be overly optimistic. We consider three aspects 
of the LHC and Oiz/3/3-decay physics not included in 
Refs. [221 [SU]: (a) the impact of SM and detector back¬ 
grounds on the significance of an LHC LNV signal; (b) 
running of the corresponding LNV effective operators 
from the TeV scale to the low-energy scale relevant to 
0://3/3-decay; and (c) long-distance contributions to the 
0:z/3/3-decay nuclear matrix element (NME). The impacts 
of these considerations are, respectively, to (a) degrade 
the significance of the LHC LNV signal for a given choice 
of LNV model parameters; (b) reduce the strength of the 
Ozz/3/3-decay amplitude relative to the inferred value of 
parameters at the high scale; and (c) enhance the NME. 
We then find that for a limited range of heavy parti¬ 
cle masses, existing Ojz/3/3-decay searches and Run H of 
the LHC may have comparable sensitivities to TeV scale 
LNV, depending on the values of the Ozz,3/3-decay nuclear 
and hadronic matrix elements. Accumulation of addi¬ 
tional data with the high-luminosity phase of the LHC 
would be necessary to achieve a reach comparable to the 
tonne-scale Ozz,3/3-decay searches. 

To be concrete, we focus on one of the simplified 
models yielding the greatest LHC reach according to 
Refs. [221 1^ - The model includes a scalar doublet S 
transforming as (1,2,1) under SU(3)cxSU(2)lxU(1)v 
and a Majorana fermion F that transforms as a SM gauge 
singlet. The interaction Lagrangian is 

/Ilnv = giQy^S + g 2 e^^UFS* + h.c. , ( 1 ) 

where L and Q are first generation left-handed lepton 
and quark doublets, respectively; d is the right-handed 
down quark; and Roman and Greek indices correspond 
to SU(2 )l and SU(3)c components, respectively. In high 
energy proton-proton collisions, the interaction Q will 
generate a final state with a same sign (SS) di-electron 
pair along with two high-pr jets. When either the S or 
F appears as an s-channel resonance, the corresponding 
cross section will be enhanced. For the low-energy Ozz/3/3- 
decay process, one may integrate out the heavy degrees 


of freedom, yielding the dimension-nine LNV interaction: 

>C!,nv = §(^i+h.c. , Oi=QT+dQT+dLL^ , (2) 

where is the lepton doublet charge conjugate field. 
Cl = glg 2 and A® = MgMp. 

We have implemented the model Q in Madgraph 
and generated events with Madevent [3T] for pp colli¬ 
sions at 14 TeV, carrying out showering, jet matching, 
and hadronization with Pythia [32] and detector simu¬ 
lation with PCS. The dominant backgrounds involve (a) 
“charge flip”, wherein one lepton from a SM opposite sign 
(OS) di-electron pair transfers most of its pp to an elec¬ 
tron of the opposite sign through conversion and (b) a 
high-pj’ jet is registered as an electron in the electromag¬ 
netic calorimeter (“jet fake”). The largest contributors 
to the charge flip background are SM production of a Z 
and virtual 7 plus jets, followed by ti production wherein 
the 5-quarks from the top decays are not tagged. For the 
jet fake background, SM multi-jet production is by far 
the leading contributor. Subdominant backgrounds in¬ 
clude diboson (WW, WZ, ZZ) plus jets. The charge flip 
background from the various aforementioned sources was 
derived by binning events in pseudo-rapidity (rj) and ap¬ 
plying the Ty-dependent charge-flip probabilities as mea¬ 
sured by ATLAS [53]. For the jet-fake background, we 
applied a medium jet-fake probability of 2 x 10 “"^ [531 IMj 
times a combinatoric factor associated with the number 
of jet-fakes in an event with JV jets. 

After imposing a set of basic selection cuts {pr. > 
20 GeV, \r]j\ < 2.8, \r]e±\ < 2.5) we find that addi¬ 

tional cuts on Ht{ iets), the scalar sum of all jet pr, 
the dilepton invariant mass, and missing energy are 
highly effective in reducing the background while main¬ 
taining the signal. A set of cuts that optimizes the signif¬ 
icance S/yJS + B is given in Table|^ The signal indicated 
is generated for Ms = Mp = 1 TeV and gi = g 2 = 0.176, 
corresponding to a Ozz/3/3-decay rate consistent with the 
present GERDA upper bound (see below). 

In order to translate the sensitivity to the parameters 
that enter the high energy process to the Ozz/3/3-decay 
rate, we evolve the operator Oi to the GeV scale using 
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the renormalization group. Between the scale /r = A 
and the weak scale v = 246 GeV, one must in princi¬ 
ple include both QCD and electroweak corrections. As 
the latter are generally considerably smaller than the for¬ 
mer, we include only the QCD corrections and continue 
the running from v to fi = 1 GeV. Under this evolu¬ 
tion, Oi will mix with three additional operators: O 2 = 
Qafji,uT~^dQa^'^T^dLL^ , O3 = QT^r^dQT^r^dLL^ , 
and O 4 = QadQ(7 >^'^T^T~^dLL^, where A = 
1 , • • - 8 denote the SU(3)c generators in the fundamen¬ 
tal representation and is the isospin raising operator. 
The corresponding anomalous dimension matrix is 


T 

7 = 


2-k 


f -8 0 0 -32/3 \ 

0 -8/3 2/9 0 

0 -48 1 -20 

V -1 0 -5/12 -19/3 y 


(3) 


The Wilson coefficients C'^ = (Ci, • • •, C 4 ) then evolve 
according to dC'/dln/r = Under this evolution, we 

find, for example, that if only Ci(/r = A) is non-vanishing 
at the high scale, then the magnitude of the Wilson co¬ 
efficients Cj(/r = 1 GeV) are: Ci = 0.203C'i(A), C 2 = 
-0.007C'i(A), C 3 = 0.266Ci(A), and C 4 = -0.055Ci(A). 

For fi below ^ 1 GeV, use of quark degrees of free¬ 
dom is no longer appropriate, so one must match the 
operators Oj onto operators built from hadronic de¬ 
grees of freedom. To that end, we follow Ref. [3S] and 
exploit the transformation properties of the Oj under 
SU(2 )l xSU(2)fl chiral symmetry. It is convenient to 
Fierz transform © 3,4 to forms in which all quark blin- 
ears are color singlets, leading to an effective coefficient 
of Oi given by 


CeS 


Ci(l GeV) - ^C 3 (l GeV) = 0 . 092(71 (A) (4) 


where we have omitted the negligible contributions from 
^ 2 , 4(1 GeV). Using the notation of Ref. [3S] we note that 
the part of Oi relevant to the decay process is 

= §1 dLe% + h.c. , (5) 

where = [eh)^ and 

02± = qRT°-qLqRT^qL ± qLT°-qRqLT^qR ( 6 ) 

with q£ = {u,d)j^ ]j. Since O 2 - is parity-odd and 
the 0^'/3/1-decay processes of experimental interest involve 
0+ —>■ 0+ transitions, we retain only the O 2 + part of 
At the hadronic level, matches onto the two 

pion-two electron operator 

-k h.c. ^ TT-TT-eLe% + h.c., 

(7) 

where = 92.2 ± 0.2 MeV is the pion decay constant 
[55] and Ah is a mass scale associated with hadronic 


matrix elements of the four quark operator Using 

the vacuum saturation and factorization approximation, 
we estimate the latter to be Ah = m‘^/{mu+md) « 2.74 
GeV for m^+ =139 MeV and -I- rud = 7 MeV [57] . 

The effective pion-electron interaction in Eq. 0 leads 
to a long-range contribution to the Qvjdjd amplitude[55]. 
Following Ref. |5S| we then obtain the following result 
for the decay rate: 



with 9c being the Gabibbo angle, I{Epis) the electron 
phase space integral 

l>E00—me 

I dEiF{Z-\-2^ Ei)F[Z-\-2^ E 2 )piEip 2 E 2 , (9) 

E 2 = Epp — El, and E(Z + 2 ,Ei^ 2 ) being factors that 
account for distortion of the electron wave functions in 
the field of the final state nucleus. The NME is given by 

= + ^ 27 ^ 7 ] r+T+ 14 /,) ( 10 ) 

where Tij = 3cji ■ pijaj ■ pij -di-dj, R = pij is the 

separation between nucleons i and j, and the functions 
Pi. 2 {\pij\) are given in Ref. [35|. Note that we have nor¬ 
malized the rate to the conventionally-used factor Gq^, 
that contains quantities associated with the SM weak in¬ 
teraction, even though the LNV mechanism here involves 
no SM gauge bosons. The rate (|^ is similarly insensitive 
to the nucleon axial vector coupling pA and the debate 
over its “quenching” in nuclei |55M5] . 

Values for Mq have been computed using the quasipar¬ 
ticle random phase approximation (QRPA) in Ref. jlSj 
for a variety of isotopes. For illustrative purposes, we 
consider the Or f3/3-deca.y of ^®Ge, for which the authors 
of Ref. [IS] give Mq = —1.99. We emphasize, however, 
that both the hadronic matching scale Ah and the NME 
Mq are presently subject to considerable theoretical un¬ 
certainties. In the case of Oi//3/3-decay mediated by the 
exchange of light Majorana neutrinos, for example, NME 
computations obtained using the nuclear shell model are 
typically a factor of two smaller than those obtained 
using QRPA. In order to illustrate the impact of both 
sources of uncertainty, we show results for two different 
values the product MoA’jj that differ by a factor of two. 

To illustrate the present and prospective reach of 
0i/f3f3-deca.y and LHG searches, we first show in Fig. 
the significance of a possible LHC observation, assum¬ 
ing Ci/A^ has the maximum value consistent with the 
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FIG. 1: Significance of a LHC e“e“ + di-jet signal as a func¬ 
tion of integrated luminosity assuming the maximum Ci/A® 
consistent with the GERDA QvjSP half-life limit. Upper and 
lower curves correspond to values of the NME Mq = —1.0 and 
— 1.99, respectively. 



FIG. 2: Present and future reach of OvPP and LHC searches 
for the TeV LNV interaction (IH as functions of the effective 
coupling ^eff and mass scale A'^(see text). Present GERDA 
exclusion and future tonne-scale sensitivity are indi¬ 

cated by upper and lower shaded regions, respectively. Darker 
shaded bands indicate impact of varying MoA'jj by a factor of 
two. LHC exclusion reach for representative integrated lumi¬ 
nosities are indicated by the solid, dashed, and dotted lines. 


present GERDA limit for ^®Ge (T 1/2 <3x10^® yr) as 
implied by Eq. Q. We see that non-observation with 
~ 735 fb~^ (~ 70 fb“^) would imply exclusion at a level 
consistent with the present GERDA limit assuming the 
larger (smaller) value of MoA'jj. The corresponding re¬ 
quirement for discovery S/y/S -|- B > 5 is ;$4.6 ab“^ (Jl 
435 fb“^). It is striking that a factor of two difference in 
MoA|j, when translated into an upper bound on Ci/A^, 
implies an order of magnitude difference in the luminos¬ 
ity needed for LHG exclusion or discovery. The ex¬ 
clusion and discovery reaches for both the LHG and a 



FIG. 3: Same as Fig. but giving LHC discovery reach. 


future, one-ton Oj^/3/3-decay as functions of A and an ef¬ 
fective coupling ^eff = C'i(A)^/^ are shown Figs. and 
1^ respectively. We use a prospective ^®Ge sensitivity of 
Ti /2 = 6 X 10^^ yr[3ni- We also show the present GERDA 
exclusion for reference. The darker shaded bands at the 
lower edges of each Oj^/3/3-decay exclusion and future sen¬ 
sitivity regions indicate the impact of varying MqA^ by 
a factor of two. From Fig. [^we observe that with > 100 
fb“^ the LHG would begin to extend the present GERDA 
exclusion for A in the vicinity of 1.4 TeV for the larger 
value of |A/o|A|^ and for a broader range of masses as¬ 
suming the smaller value. As indicated by Fig. the 
opportunities for discovery with 300 fb“^ appear more 
limited, even under the assumption of the smaller nu¬ 
clear and hadronic matrix elements. However, the high 
luminosity phase of the LHG with 3 ab“^ could open 
the possibility for discovery over a range of masses that 
depends on the value of A/oA|^. 

From the standpoint of the LHG, this conclusion is not 
as optimistic as obtained in Refs. [Miisn], as the reach 
of the tonne-scale Oz^/3/3-decay experiments appears to 
exceed that of the high-luminosity LHG over nearly the 
entire range of parameter space considered. It is, nev¬ 
ertheless, interesting to compare the prospects for both 
Oj^/3/3-decay and the LHG, as observation of a signal in 
both experiments is possible and would point to the exis¬ 
tence of TeV scale LNV interactions. Reducing the /3,0- 
decay nuclear and hadronic matrix element uncertainties, 
as well as refining the estimates of jet-fake and charge flip 
backgrounds at the LHG, would clearly clearly sharpen 
the implications of this comparison. 
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